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Correlations between the behavior of the nuclear symmetry energy, the neutron skins, and the percent- 
age of energy-weighted sum rule (EWSR) exhausted by the Pygmy Dipole Resonance (PDR) in ^*Ni 
and ^^^Sn have been investigated by using different Random Phase Approximation (RPA) models for 
the dipole response, based on a representative set of Skyrme effective forces plus meson-exchange effec- 
tive Lagrangians. A comparison with the experimental data has allowed us to constrain the value of the 
derivative of the symmetry energy at saturation. The neutron skin radius is deduced under this constraint. 
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One of the interesting problems presently receiving par- 
ticular attention, is that of the size of the neutron root- 
mean-square (r.m.s.) radius in neutron-rich nuclei. In fact, 
this quantity is related to the isospin-dependent part of the 
nuclear equation of state (EOS) which, in turn, has rele- 
vant implications for the description of neutron stars. At 
present, there is an enormous effort aimed at determining 
the parameters that govern the asymmetric matter EOS, us- 
ing both experimental and theoretical tools. Review papers 
have been devoted to this topic 111121. 

The energy per particle in a nuclear system characterized 
by a total density p (sum of the neutron and proton densities 
Pn and pp), and by a local asymmetry 6 = {p„ — Pp) / p, 
is usually written as 

^ip,5) = ^{p,5 = 0) + Sip)6'. (1) 

Odd powers of 5 are forbidden by the isospin symmetry 
and the term proportional to 6'^ is found to be negligible. 
The above equation defines the so-called symmetry energy 
S{p). Determining values of the symmetry energy at var- 
ious densities of interest for nuclear structure, nuclear re- 
actions, and astrophysics, is one of the great challenges for 
the physics community. 

Information on the symmetry energy can be obtained 
from various sources, none of them being so far conclu- 
sive by itself. A direct correlation between the neutron skin 
thickness AR and the derivative of the symmetry energy at 
saturation, has been found in Refs. |[3j|4l. The derivative of 
the symmetry energy at saturation is related to the widely 
used "slope" parameter L by 

S'(.p)l=po = ^- (2) 

The symmetry energy at saturation, S{po), is denoted by 
a4 or J: we shall use the symbol J in what follows. No 



measurement of the neutron skin is available which is accu- 
rate enough to constrain the slope parameter L. The prop- 
erties of the isovector Giant Dipole Resonance (IVGDR) 
[5 1, of the low-lying electric dipole excitation (the so-called 
Pygmy Dipole Resonance, PDR) [61, and of the charge- 
exchange spin-dipole strength [7| have been suggested as 
constraints. In addition, by means of heavy-ion collisions 
the symmetry energy has also been probed at subsaturation 
densities (0.4 < p/pa< 1.2). In Ref. |[8l, isospin diffusion 
data from the collision between ^^^Sn and ^^''Sn have been 
analyzed using a transport model in which the momentum- 
dependent symmetry potential enters as one of the main 
ingredients. The same data, together with the double ratios 
of neutron and proton energy spectra, have been analyzed 
within a different kind of transport model in Ref. [9|. We 
should also mention that another analysis of isoscaling data 
has been reported in Ref. fflOl . Finally, in the work re- 
ported in Ref. iflTI a range of values for L is inferred from 
the analysis of data of radii from antiprotonic atoms. 

One of the motivations of the present work lies in the 
consideration that the values of L extracted from the PDR 
in "^sn (between ^ 30 and 60 MeV) [W] and from the 
GDR in ^'^^Pb [12|, are smaller than those deduced from 
the analysis of the heavy-ion collisions. More precisely, 
the values of L deduced with the two different approaches 
overlap only in a small interval (cf., e.g.. Fig. 3 of [9|). We 
would like to pursue in this paper an analysis which is more 
general than the one performed in Ref. [6|, by considering 
PDRs in two different mass regions and a variety of the- 
oretical models, both nonrelativistic and relativistic. Both 
classes of mean-field models are successful in describing 
the nuclear ground states and of many of the excited states 
(for a review, see Ref. ifTSll ). Our goal is to see whether 
consistency among different ways of extracting the slope 
parameter L can be achieved: as a result, one should also 
expect to be able to better pin down the values of the neu- 
tron skin radii. 
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Progress in the study of the density dependence of the 
symmetry energy and the neutron radii through the PDR 
requires more work in two directions. The first is to have 
more data, particularly for unstable neutron rich-nuclei 
characterized by a sizeable dipole strength in the low- 
energy region. The second is a more comprehensive the- 
oretical analysis of the PDRs based on different calcula- 
tions both of nonrelativistic and relativistic RPA types. In 
this paper this is realized by taking advantage of a recent 
experimental datum on the PDR in the neutron-rich ^^Ni 
nucleus Id 4-1 . In particular, the present analysis has aimed 
(i) at finding evidence of a correlation between the values 
of L and the EWSR exhausted by the PDRs when they are 
calculated using many Skyrme parameter sets and covari- 
ant effective Lagrangians, (ii) at inferring values for the 
neutron skin radii of ^^Ni, ^^^Sn and ^°^Pb, and finally (iii) 
at comparing our deduced value of L with the other ones 
existing in the recent literature. 

The first step consisted in finding correlations between 
the properties of the PDR and the symmetry energy. Ran- 
dom Phase Approximation (RPA) calculations of the dipole 
strength have been carried out. Our implementation based 
on nonrelativistic Skyrme forces is fully self-consistent and 
discussed, e.g., in Ref. ifTSl . The Hartree-Fock (HF) equa- 
tions are solved in a radial mesh extending up to ~4 times 
the nuclear radius. The continuum is discretized by using 
box boundary conditions. The model space is large enough 
so that the well-known double commutator EWSR is ex- 
hausted by at least 96%. We have employed 19 different 
Skyrme sets which can be said to constitute a quite repre- 
sentative ensemble. All of them have an associated value of 
the nuclear incompressibility K^a lying in the interval 210- 
270 MeV [16|. We do not provide here the original refer- 
ences in which the parameter sets have been introduced: 
they can be found in jJI \Vf\. We have checked that our 
results do not change appreciably if we take out, or add, a 
few Skyrme parameter sets to our ensemble. The relativis- 
tic calculations are based on the well-known relativistic 
mean field (RMF) theory plus the self-consistent relativis- 
tic RPA (RRPA) as described in Refs. USIUl. We have 
employed 7 different parametrizations for the non-linear, 
meson-exchange effective Lagrangian. In the calculations, 
box boundary conditions are used: the box radius is set at 
30 fm and the radial mesh is 0. 1 fm. The model spaces for 
particle-hole and antiparticle-hole are cut at energy Scutoff 
= 1039 MeV and -939 MeV, respectively. The references 
for the seven parameter sets can be found in ||20l|2]l. 

We have found a rather good correlation between the pa- 
rameter L and the percentage of EWSR associated with the 
PDR. In the theoretical calculations, we consider the whole 
part of the low-energy region where the strength is not neg- 
ligible. In the nucleus ^^Ni the PDR is associated, as a 
rule, with a well-defined peak between 9 and 1 1 MeV, to be 
compared with the experimental finding of [141, that is, 1 1 
MeV. In few cases the strength is more fragmented and/or 
at lower energy. We display two typical dipole strength 



distributions in Fig. [T] the separation between PDR and 
GDR regions looks quite clear. In the nucleus ^'^^Sn, the 
peak of the PDR is between 7.5 and 9.5 MeV. The exper- 
imental peak energy is 9.8 MeV ||6l. The percentages of 
EWSR are defined in this work with respect to the classical 
Thomas-Reihe-Kuhn (TRK) value, and vary between 1% 
and 10%. In general, the relativistic Lagrangians provide 
larger values for this latter quantity. The PDR energies they 
provide in the case of ^^^Sn are also about 1 MeV lower 
than the experimental value: therefore, trying to constrain 
the symmetry energy by using the correlation between the 
PDR energy and the value of S at p=0.l fm~^ (plus the 
experimental datum) was attempted in ll20l . but it was only 
possible by means of extrapolation. 

In the upper part of Fig. |2] the correlation between the 
percentage of the EWSR and L is shown for both nuclei 
^*Ni and ^^■^Sn. The straight lines correspond to linear fits. 
We have considered the measured values of the EWSR per- 
centage, and deduced a range of acceptable values for L, 
by taking care both of the experimental error and of the 
error associated with the fit (the latter being almost neg- 
ligible with respect to the former). Our results are more 
general than those presented in Ref. [6 1 since we consider 
two different nuclei and many different mean field models. 
Although we do not include in our analysis all classes of 
mean field models, we try nonetheless to avoid, as much 
as we can, possible sources of bias since we avoid restrict- 
ing to Skyrme sets fitted by the same group with the same 
protocol. In fact, our sets span a broad range of possible 
values associated with nuclear matter quantities. 

In the case of ^^Ni the measured value of the EWSR per- 
centage is 5%ibl.5%. The error includes the uncertainty 
related to the quantities used to deduce the number from 
the measurement. It should be noted that the dominating 
uncertainty (still within 30% of the average value) is that 
related to the choice of the level density value entering the 
evaluation of the branching for gamma emission. We have 
used different level densities obtained by means of either a 
Shell Model Monte Carlo (SMMC) calculation for this nu- 
cleus 11221 . or global Hartree-Fock-Bogoliubov (HFB) cal- 
culations ll23l l24ll : the largest span goes from 3.5% ob- 
tained using ll23l to 6.5% using ll24ll . 

Our result is that the slope parameter L is constrained 
to be in the interval 50.3-89.4 MeV or 29.0-82.0 MeV, if 
we use either the ^^Ni results, or the ^"^^Sn results (cf. the 
lower left panel of Fig. |2|. The weighted average, L = 
64.8ibl5.7 MeV, is displayed in the lower right panel of 
Fig. [2](it corresponds to the shaded box). In this panel, the 
correlation of J and L is provided, so that we can deduce 
our best value of J which is 32.3ibl.3 MeV. This value is in 
very good agreement with the value 32.0ibl.8 MeV which 
is reported in Ref. ||6j|. The parametrizations of S{p) found 
in Refs. lead to J = 31.6 MeV. Moreover, our re- 

sult for J overlaps well with the ranges obtained in Refs. 
m (30.2-33.8 MeV) and [25| (31.5-33.5 MeV) (cf. also 
ll26l ). From the theoretical point of view, we can consider 
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very satisfactory that our result for L coincides almost ex- 
actly with the value of 66.5 MeV extracted from Bruckner- 
Hartree-Fock (BHF) calculations in uniform matter that 
employ realistic two-body and three-body forces f27l . 

The next step is to use the L value obtained from the 
PDR computed data points in ^*Ni and ^'^^Sn in order to 
deduce the neutron skin thickness Aii. First, one can 
note that the correlation between L and Ai?, when the 
two quantities are calculated using the models already de- 
scribed, is quite good (cf. Fig. [3]l. If one imposes the 
value of L to be in the interval 64.8ibl5.7 MeV, one ob- 
tains for the skin thickness Ai?=0.200±0.015 fm for ^^Ni, 
Ai2=0.258±0.024 fm for "^Sn^ and Ai?=0.194±0.024 
fm for ^"'^Pb. These numbers are stable if one tries to con- 
strain them by using the L value from ^'^Ni only, or ^■^^Sn 
only, instead of the weighted average. It should also be 
noted that the values associated with IS.R, both for ^'^^Sn 
and ^°^Pb, are in good agreement with the results reported 
in Ref. [6|. This gives us further confidence on the value 
of the neutron skin of ^^Ni which is determined for the first 
time through the present analysis. We should recall that 
the possibility to extract Aii directly from measurements 
of the spin-dipole strength has been discussed |28|. For a 
thorough discussion of the exensive literature appeared in 
previous decades on this subject we refer the reader to Ref. 

In Fig. |4] we show the comparison of the values of L 
found in our analysis with those found with other analy- 
sis and/or other methods. The main point is that the result 
for L extracted from the PDR in ^■^^Sn is compatible with 
the one from Ref. [6|; however, combining in our analy- 
sis the two PDRs of both ^''^Sn and ^^Ni, we are able to 
shift the range of L to larger values and to reduce the un- 
certainty. This solves, to a good extent, the problem that 
the result from Ref. [6| was not overlapping significantly 
with the results obtained by the different analysis of heavy- 
ion collisions. In the lower panel of Fig. |4j one can see 
that our present finding has a remarkable overlap with the 
results of most of the other proposed methods to extract 
L, that involve not only different methodologies but also 
very different observables. We can conclude that our more 
general analysis of the extraction of the slope parameter L 
from the PDR is able to provide a firm result. Another im- 
portant side result of our work is that we are able for the 
first time to propose a value for the neutron skin thickness 
of the neutron-rich ^^Ni isotope. More PDR data in other 
mass regions and/or in long isotopic chains are desirable to 
increase the predictive power of our procedure. This could 
lead to determine quite accurately quantities such as the 
neutron radii, and the parameters governing the density de- 
pendence of the symmetry energy, that are fundamental for 
nuclear physics and for their implications in the study of 
neutron stars. 
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FIG. 1: Two typical dipole strength functions calculated in the 
nucleus ''^Ni. A nonrelativistic and a relativistic example are 
shown in panels (a) and (b) in which, respectively, the Skyrme 
force SkI3 and the NL3 parametrization of the effective RMF 
Lagrangian have been used. The sharp RPA peaks are averaged 
by using Lorentzian functions having 1 MeV width. 
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FIG. 2: (Color online) In panels (a) and (b), the correlation be- 
tween L and the percentage of TRK sum rule exhausted by the 
PDR, respectively in ^^Ni and ^^^Sn, is displayed. The computed 
data points are labelled, here and in what follows, by numbers. 
The correspondence with the parameter sets used is: l=v090, 
2=MSk3, 3=BSkl, 4=vll0, 5=vlOO, 6=SkT6, 7=SkT9, 8=SGII, 
9=SkM*, 10=SLy4, ll=SLy5, 12=SLy230a, 13=LNS, 14=SkMP, 
15=SkRs, 16=SkGs, 17=SK255, 18=SkI3, 19=SkI2, 20=NLC, 
21=TM1, 22=PK1, 23=NL3, 24=NLBA, 25=NL3+, 26=NLE. 
The straight lines correspond to the results of the fits. In panel 
(c) we show the same straight lines displayed in (a) and (b), to- 
gether with the correlation coefficient r and the constraints from 
experiments f?", 141. In panel (d) the correlation between L and 
J is shown. The box corresponds to the value of L deduced from 
the weighted average of the two values extracted from 
"2Sn. 
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FIG. 3: (Color online) Panels (a), (b) and (c) display the cor- 
relations between the neutron skin thickness AT? and the slope 
parameter L, in the case of the three nuclei analyzed in this work. 
The convention is the same as in the previous figure. Under the 
constraint for L emerging from our analysis [shaded area in panel 
(d)], the values displayed for the neutron skin thickness in the 
three nuclei are obtained. 



7 



20 40 60 80 100 120 140 

(a) PDR analysis 

this worlc: PDR in ™Sn I 



this work^PDRin "Ni 
Klim kiewicz et a l. 2007 

_i I , I , I , \ I , , I I I , \ i_ 

(b) 

this work: different PDR,^ 



All methods 



Tsang et al. 2009 
Shetty et al. 2007 

Chen et al. 2005 
Danielewicz 2003 

Danielewicz, Lee 2009 



Cent elles et. al. 2 009 
Klimkiewicz et al. 2007 



20 40 60 80 100 120 140 
L [MeV] 



FIG. 4: (Color online) Comparison between the values of L ex- 
tracted in the present work, and those from the existing literature. 
In panel (a), we compare our separate results from the PDRs of 
^^Ni and ^^^Sn with the result of Klimkiewicz et al. \6\. In panel 
(b), we compare with values extracted from completely different 
kind of analysis: Tsang et al. |9|, Shetty et al. 1 10|, Chen et al. 
ifSi , Danielewicz |26|, Danielewicz and Lee [25 J , Centelles et al. 
mil , and Klimkiewicz et al. |[6|. 



